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ABSTRACT 

The origin and evolution of magnetic fields in the Universe is still an open question. 
Their observations in galaxies suggest strong magnetic fields already at high redshift 
as well as at present time. However, neither primordial magnetic fields nor battery 
processes can account for such high field strengths, which implies the presence of 
a dynamo process with rapid growth rates in high-redshift galaxies and subsequent 
maintenance against decay. 

We investigate the particular role played by feedback mechanisms in creating 
strong fiuid turbulence, allowing for a magnetic dynamo to emerge. Performing mag- 
netohydrodynamic simulations of isolated cooling gas halos, we compare the magnetic 
field evolution for various initial field topologies and various stellar feedback mech¬ 
anisms. We find that feedback can indeed drive strong gas turbulence and dynamo 
action. We see typical properties of Kolmogorov turbulence with a kinetic en¬ 

ergy spectrum, as well as a small-scale dynamo, with a magnetic energy spectrum 
predicted by Kazantsev dynamo theory. We also investigate simulations with a final 
quiescent phase. As turbulence decreases, the galactic fountain settles into a thin, 
rotationally supported disk. The magnetic field develops a large-scale, well-ordered 
structure with even symmetry, which is in good agreement with magnetic field obser¬ 
vations of nearby spirals. Our findings suggest that weak initial seed fields were first 
amplified by a small-scale dynamo during a violent, feedback-dominated early phase 
in the galaxy formation history, followed by a more quiescent evolution, where the 
fields have slowly decayed or were maintained via large-scale dynamo action. 

Key words: galaxies: magnetic fields - methods: numerical - MHD - turbulence 


1 INTRODUCTION 

Large-scale magnetic fields are measured with strengths of 
up to several fiG in nearby galaxies (Beck et al. 1996), and 
possibly even higher field strength have been detected in 
earlier galaxies at high redshift (Bernet et al. 2008). The 
preferred theory to explain their origin is based on the early 
generation of seed fields at the epoch of cosmic re-ionisation, 
through the microscopic process known as ’’Biermann bat¬ 
tery” (Biermann 1950; Naoz & Narayan 2013), later ampli¬ 
fied during the galaxy formation era, through the large-scale 
galactic dynamo process (Parker 1970; Brandenburg & Sub- 
ramanian 2005) . The Biermann battery is likely to generate 
fields at the level of 10“^° G, leaving to the galactic dynamo 
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process more than 14 orders of magnitude of field amplifica¬ 
tion during the 10 Gyr of cosmic evolution. The challenge for 
galactic dynamos is even more severe, if one considers that 
strong fields are already in place at high redshift (Widrow 
2002), and are probably even stronger than they are today 
(Bernet et al. 2008). 

Successful theoretical models for large-scale galactic 
dynamos report exponential growth rates of the order of 
r ~ 0.01 to O.IU, where U is the galactic disk rotation rate 
(Pariev et al. 2007). For typical, present day spirals, this 
translates into e-folding amplification time scale of roughly 
1 Gyr, making the task of amplifying the field over 14 or¬ 
ders of magnitude virtually impossible. One noticeable ex¬ 
ception is the cosmic-ray-driven dynamo proposed by Parker 
(1992) and simulated by Hanasz et al. (2004), leading to a 
measured growth rate F Cl, although the numerical ex- 
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periment was performed over only a relatively limited time, 
since the reported magnetic field amplification was over only 
3 orders of magnitude (Hanasz et ah 2004). On the theoret¬ 
ical side, classical mean field dynamos are plagued by the 
catastrophic a-quenching effect, leading to very low satu¬ 
rated values for the large-scale magnetic field (Kulsrud & 
Anderson 1992; Vainshtein V Cattaneo 1992), owing to the 
strict conservation of magnetic helicity in a closed system. 
A possible solution to this problem is the effect of galactic 
winds, that could drag the magnetic field lines outside of 
the dynamo-active disk, therefore alleviating the aforemen¬ 
tioned quenching issue (Del Sordo et al. 2013). 

The theory of galaxy formation has significantly evolved 
over the past decade, with the ever increasing role of feed¬ 
back processes (Scannapieco et al. 2012) and their associated 
galactic winds (Oppenheimer V Dave 2006), together with 
the dominant accretion mechanism through cold streams 
(Keres et al. 2005; Ocvirk et al. 2008; Dekel et al. 2009). 
On the observational side, galactic winds are indeed ubiqui¬ 
tous in star bursting local galaxies (Martin 1999), but also in 
many “normal” high redshift galaxies (Steidel et al. 2010). 
One of the most spectacular observational constraints on 
galaxy formation theories was obtained by matching the 
stellar mass of the central galaxies to their parent halo mass 
(Behroozi et al. 2013; Moster et al. 2013). This has led theo¬ 
rists to consider much stronger feedback processes, in order 
to regulate star formation throughout cosmic time, espe¬ 
cially at high redshift, when the star formation efficiency 
was so low (Agertz et al. 2013; Hopkins et al. 2014; Roskar 
et al. 2014; Wang et al. 2015). 

In this rather violent, feedback-dominated scenario, 
dwarf galaxies play a very important role. They are the dom¬ 
inant galaxy population at high redshift, probably responsi¬ 
ble for the cosmic re-ionisation (Kimm V Cen 2014). They 
are also the progenitors of the Milky Way satellites, which 
are useful laboratories to test our current galaxy formation 
paradigm. For the latter, violent feedback mechanisms have 
also been invoked to explain the absence of cusp in the dark 
matter density profile, and the presence of a dark matter 
core in low surface brightness galaxies (de Blok et al. 2001). 
Cosmological simulations of dwarf galaxies have been per¬ 
formed with strong feedback recipes, conhrming in this case 
the formation of a dark matter core (Governato et al. 2010, 
2012). Recently, we have also performed idealised simula¬ 
tions of an isolated, cooling gaseous dwarf halo, obtaining, 
in this well-controlled numerical experiment, the formation 
of a dark matter core (Teyssier et al. 2013). The dark matter 
core formation mechanism is now well understood (Pontzen 
V Governato 2012). It is due to repeated, energetic feed¬ 
back events due to many supernovae explosions, leading to 
violent oscillations of the gravitational potential, due to the 
large gas mass variations within the central kilo parsec of 
the galaxy. A possible observational signature of this effects 
is a typical, oscillatory star formation history, mimicking a 
breathing mode in the gas distribution (Kauffmann 2014). 

In this paper, we want to study the impact of a strong 
feedback scenario on the growth of magnetic fields in dwarf, 
as well as in larger galaxies. The velocity field on both 
small and large scales, resulting from repeated giant feed¬ 
back events, can have a direct influence on the growth of the 
magnetic energy. Indeed, supernovae explosions in the Milky 
Way have been considered for quite a long time as a source 


of helical gas motions, promoting a large-scale a-dynamo 
in the Galaxy (Ferriere 1992). The Milky Way is however 
a rather quiescent galaxy, with moderate supernovae activ¬ 
ity. In this paper, we are considering feedback-dominated 
galaxies, with high star formation rates and violent turbu¬ 
lent motions, together with large-scale galactic fountains or 
winds. 

Several simulations including magnetic fields have been 
performed recently in the context of galaxy formation (Wang 
& Abel 2009; Dubois & Teyssier 2010). These simulations, 
based on the popular “cooling halo” numerical set-up, have 
achieved only moderate magnetic field amplification. The 
important property of these simulations is the absence of 
feedback (Wang & Abel 2009), or the relative weakness of 
the feedback recipe used at that time (Dubois Sz Teyssier 
2010 ). 

A first exception is the simulation reported in Beck et al. 
(2012), based on a MHD version of the SPH code GAD¬ 
GET with divergence cleaning. They observed a fast expo¬ 
nential growth of the magnetic field, which they attribute 
to a small-scale dynamo. Feedback processes were included 
through an effective Equation-of-State (EoS), without any 
explicit source of turbulence in these relatively smooth, 
thermally-supported flows. These authors however reported 
very strong growth rates, with e-folding times as small as 
10 Myr, although analytical estimates based on small-scale 
dynamo theory predicted e-folding times closer to 100 Myr. 

A second exception is the recent simulation reported 
in Pakmor & Springel (2013), using the new Magneto- 
Hydrodynamics (MHD) solver developed for the AREPO 
code (Pakmor et al. 2011), where strong magnetic field am¬ 
plification has also been observed, although, here again, stel¬ 
lar feedback effects were not considered explicitly, but only 
indirectly as a modified thermal EoS, leading to the forma¬ 
tion of relatively smooth, two-dimensional flows, in which 
dynamo amplification is in principle notoriously difficult to 
obtain. 

In the present paper, we will use a similar set-up as in 
all those previous studies, namely a cooling isolated gaseous 
halo, considering simulations with (but also without) strong 
stellar feedback. We will use the Adaptive Mesh Refinement 
code RAMSES (Teyssier 2002), adopting the “Gonstrained 
Transport”, strictly divergence-free-preserving, MHD solver 
presented in Teyssier et al. (2006) and in Eromang et al. 
(2006). The paper is organised as follows: In section 2, we 
will present our numerical methods, both in terms of galaxy 
formation physics and magnetic field modelling. In section 3, 
we describe our initial conditions for the isolated, magne¬ 
tised cooling halo. In section 4, we present our main results, 
outlining the difference between the feedback and the no¬ 
feedback cases. Einally, in section 5, we discuss the implica¬ 
tions of our results in the context of galactic dynamo theory, 
as well as possible further studies to confirm and broaden 
our findings. 


2 NUMERICAL METHODS 

We have performed MHD simulations of isolated, cooling 
haloes, using the Adaptive Mesh Refinement (AMR) code 
RAMSES (Teyssier 2002). These simulations feature a colli¬ 
sionless fluid (for dark matter and stars) and a magnetised 
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gaseous component, coupled through gravity. In this sec¬ 
tion, we describe the simulation technique used to follow 
the evolution of our isolated halo. First, we describe in de¬ 
tails our AMR implementation for solving the ideal MHD 
equations, together with simple test cases to show that it 
works as intended in the context of galactic dynamo. We 
then describe the adopted galaxy formation physics, such 
as gas cooling, metal enrichment, star formation and stellar 
feedback, leading to what we believe to be a realistic model 
of the interstellar medium (ISM). 

2.1 Ideal MHD solver 

The equations that we solve are the ideal MHD equations 
(written here without gravity and cooling source terms for 
the sake of simplicity) 


dtp + V • {pu) 

= 0 

( 1 ) 

dt{pu) + V • {puu^ - BB^ + Ptot) 

= 0 

( 2 ) 

dtEEV- [{E + Ptot)u - {u • B)B] 

= 0 

( 3 ) 

dtB — V X (it X H) 

= 0 

( 4 ) 


where p is the gas density, pu is the momentum, B is the 
magnetic field, E — \pu^ + £ + is the total energy, 
and £ is the internal energy. The total pressure is given by 
Plot — P -\- where we assume a perfect gas equation of 
state P = (7 — l)e. The system of equations is completed 
by the soloinoidal constraint 

V • H = 0. (5) 

The code is grid-based with a tree-based adaptively re¬ 
fined mesh. The equations are solved using the second-order 
unsplit Godunov scheme based on the MUSCL-Hancock 
method. We chose the HLLD Riemann solver with the Min- 
Mod slope limiter for the hydro variables which are cell- 
centred. The magnetic field on the other hand is treated 
as a face-centered variable. This allows the use of the Con¬ 
strained Transport (CT) method to advance the induction 
equation (Equation 4) in time, which preserves the diver¬ 
gence of the magnetic field to the numerical precision level 
(Teyssier et al. 2006). The CT method involves a spatial in¬ 
terpolation of the EMF on the cell edges for the predictor 
step and solving a 2D Riemann problem for the corrector 
step. Eor the 2D problem, we use the HLLD solver as well 
and for the magnetic field in general, the MonCen slope lim¬ 
iter. 

Boundary conditions were chosen to allow for free out¬ 
flow. For the 5 hydro variables, this is done by imposing a 
vanishing gradient at the domain boundary (zero-gradient 
method). The same can be applied to the transverse mag¬ 
netic field component parallel to the boundary face Py, but 
would cause a non-zero divergence of the magnetic field if 
applied to the normal component B±_ which is perpendic¬ 
ular to the face. Instead, we use a linear interpolation for 

so that V ‘ B — 0. Note that this method can cause 
an inward Poynting flux which transports magnetic energy 
from the outside into the computational domain. Since the 
magnetic field at the border is many orders of magnitude 
weaker than the average, this does not contribute signifi¬ 
cantly to the overall magnetic energy evolution (see Dubois 
& Teyssier 2010). 

Special care needs to be taken also when refining and 


de-refining cells, in order to enforce the V-H = 0 constraint, 
when interpolating the magnetic field. A solution to this 
problem within the CT framework has been proposed by 
Balsara (2001) and Toth & Roe (2002), and we adopt it 
here for newly refined cells, but also for temporary ghost 
cells used to set proper boundary conditions at coarse-fine 
level boundaries. 

In the context of galactic dynamos, it is worth mention¬ 
ing that our code has been tested extensively against well- 
known flows triggering fast dynamos, such as in the ABC 
flow (Galloway & Frisch 1986; Childress & Gilbert 1995) or 
in the Ponomarenko dynamo (Ponomarenko 1973). We have 
shown in Teyssier et al. (2006) that our numerical scheme for 
the ideal MHD equations was in fact slightly resistive, with, 
for a regular Cartesian grid, a numerical magnetic Reynolds 
number roughly inversely proportional to the square of the 
number of grid points per box length. This scaling is to be 
expected for second-order schemes and smooth flows. In the 
context of AMR and highly complex, turbulent flows, de¬ 
termining the exact effective numerical Reynolds number of 
the simulated flow is impossible. Qualitatively, though, it is 
important to bear in mind that magnetic reconnection and 
other diffusive processes occur in the simulation at a typical 
scale probably very close to the grid scale. This scale plays a 
very important role in dissipating the kinetic energy of the 
turbulence, and also controls the magnetic energy losses due 
to reconnection or (numerical) Ohmic dissipation. 

2.2 Cooling and star formation 

In addition to solving the ideal, self-gravitating MHD equa¬ 
tions, we also include many physical processes relevant to 
galaxy formation. One of the key physical ingredient is gas 
cooling, which leads the hot, initially hydrostatic halo gas 
to loose pressure support and to condense in the centre as a 
centrifugally supported disc. When this atomic gas of 10^ K 
is allowed to cool even more due to fine-structure metal line 
cooling or molecular cooling, the disc fragments into dense 
clumps, leading to the formation of a turbulent, multiphase 
medium. To model gas cooling, we use standard H and He 
cooling processes with an additional contribution from met¬ 
als based on Sutherland & Dopita (1993) for temperatures 
above 10^ K and metal fine-structure cooling below 10^ K, 
as in Rosen et al. (1995). The metallicity, denoted as Z, is 
modelled as a passive scalar, representing the mass fraction 
of atoms heavier than Helium in the gas. It is initialised to 
Zini = O.O 5 Z 0 in the halo, mimicking molecular Hydrogen 
cooling in a zero metallicity gas. The metallicity is increased 
further by supernova feedback events, using a metal yield of 
2 /= 0 . 1 . 

Our refinement strategy is based on a quasi-Lagrangian 
approach: each cell is refined if it contains more than 8 
dark matter particles or if its baryonic mass (including 
gas and star particle mass) exceeds 8 x mres, where mres 
is the adopted mass resolution of the simulation. Refine¬ 
ment are performed recursively, on a cell-by-cell basis, until 
the adopted maximum level of refinement is reached (noted 
^max). It is crucial for astrophysical simulations to resolve 
spatially the Jeans length (Truelove et al. 1997). Requiring 
that the Jeans mass is resolved by at least 64 mass resolution 
elements, Mj = 64mres, and adopting a realistic minimum 
temperature for the gas, noted Tj, one can compute the cor- 
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responding Jeans length, and require it to be resolved by 4 
cells, Aj = 4Axmin- We can then determine the maximum 
required level of refinement corresponding to the adopted 
mass resolution TTires- To prevent the gas from accumulating 
and locally violating the Jeans length criterion, we also add 
an artificial pressure floor, 

Pj = (4Aa;^in)"Tp2 (6) 

Try 

so that the gas density will never significantly exceed a typ¬ 
ical value nj given by the relation feTj = Tj(nj)/nj. 

Stars are treated as collisionless particles which are cre¬ 
ated stochastically from the gas according to a Schmidt law 
(as in Rasera & Teyssier 2006) 



if the local density pgas is above a threshold density p* = 
n^cmH. In this paper, we always choose the star formation 
threshold density to be equal to the previously defined Jeans 
density nj. The star formation efficiency per free-fall time is 
always set to e* = 0.01; this value is based on observations 
of nearby molecular clouds (Krumholz & Tan 2007). Cre¬ 
ation of stellar particles is a local random Poisson process 
with Poisson parameter A = where At is the 

simulation time-step and 

771* = 77* (AXmin) (8) 

the mass of the resulting stellar particle, which is equal to 
the smallest cell mass at the density threshold. For each sim¬ 
ulation, precise numbers for all these parameters are given 
in Table 1. 


2.3 Stellar feedback 

In this paper, we explore the consequences of strong feed¬ 
back scenarios on the amplification of the magnetic field in 
dwarf and Milky-Way-sized galaxies. The proper modelling 
of stellar feedback mechanisms, such as supernovae explo¬ 
sions, photo-ionised bubbles or infrared radiation in dusty 
environment is subject to intensive research throughout the 
ISM and galaxy formation literature. Understanding in de¬ 
tails these various processes goes far beyond the scope of this 
paper. Our goal is merely to use various phenomenological 
recipes to model such feedback mechanisms very crudely, 
and produce dynamical properties that we believe are rele¬ 
vant for high-redshift galaxies, the most important one be¬ 
ing the gas velocity field, highly turbulent, explosive and 
fountain-like, which could result in a fast magnetic dynamo. 


2.3.1 Supernovae feedbaek 

For this purpose, we used a numerical model for supernovae 
feedback developed in the context of dwarf galaxies evo¬ 
lution, and that turned out to lead to the formation of a 
dark matter core (Teyssier et al. 2013). The main ingredient 
is the use of a non-thermal energy variable, and its associ¬ 
ated pressure, treated as a passively advected scalar quantity 
eturb = peturb, which represents various small-scale, non- 
thermal energies released by supernovae (e.g. turbulence. 


magnetic fields or cosmic rays). The evolution of this non- 
thermal energy is specified by 

= Rni - (9) 

^diss 

where the dissipation time scale is fixed to tdiss = 20 Myr 
and the energy injection per supernovae is set by 

^inj = p^TjSN • 10^° erg/M^ (10) 

where the mass fraction in massive stars is set to 77 SN = 0.1 
and the local star formation rate p* is set by our adopted 
Schmidt law. For details about the implementation, we refer 
to Teyssier et al. (2013) 

2.3.2 Radiation feedbaek 

Because supernovae are probably not energetic enough to 
trigger strong winds in Milky-Way-sized galaxies, it has been 
proposed recently to consider stellar radiation as an addi¬ 
tional feedback mechanism (Murray et al. 2010). Interstellar 
dust indeed absorbs UV photons, much of it being subse¬ 
quently re-emitted as thermal radiation in the infrared band. 
This radiation will transfer momentum to the gas through 
radiation pressure (Murray et al. 2010; Hopkins et al. 2012; 
Agertz et al. 2013; Roskar et al. 2014) In this paper, we 
consider radiation feedback only for Milky-Way-sized galaxy 
simulations. We use here again a very crude model to cap¬ 
ture the energy from the stellar UV radiation, using a simple 
escape probability model as 

Euv = Urad [1 - exp(-A7uvPdust Ax)] (11) 

with the dust mass density is assumed to be pdust = ^Pgas 
(here, Z denotes the gas metallicity). The dust opacity at 
0.1 pm is taken to be /^uv = lOOOcm^/g (Draine V Li 2007), 
and the total energy released during the first 10 Myr of a 
10 M© progenitor Urad = 10^^ erg/M©. The same cell is then 
assumed to absorb the energy in the infrared 

Em = Euv [1 — exp(—/^iRpdust Ax)] (12) 

where Km is the dust opacity in the IR band, which is a free 
parameter in our feedback implementation, usually around 
10 cm^/g (Draine V Li 2007; Semenov et al. 2003). The 
energy Em is added in the non-thermal energy equation for 
Uinj in the supernova feedback so that it contributes to eturb- 
Details about the implementation can be found in Roskar 
et al. (2014). Here again, we would like to stress that our 
goal is not to study in great details realistic feedback mecha¬ 
nisms, but rather to generate galactic velocity fields in qual¬ 
itative agreement with high-redshift galaxies and their asso¬ 
ciated strong outflows, in the context of galactic dynamos. 


3 INITIAL CONDITIONS 

We have performed a series of non-cosmological simulations 
of isolated halos in hydrostatic equilibrium, varying the ini¬ 
tial set-up with two different halo sizes (a typical dwarf and 
a typical Milky-Way) and two different initial magnetic field 
topologies (dipole and quadrupole), in addition to the var¬ 
ious options for stellar feedback that we have discussed in 
the previous section. We will now describe more precisely 
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our initial set-up, and a summary of the various run param¬ 
eters is given in Table 2. We refer the interested reader to 
Teyssier et ah (2013) for a more detailed analysis of the non- 
magnetohydrodynamic properties in the dwarf halo case. 


3.1 Initial halo 

Our initial halo follows a Navarro et ah (1997) (here¬ 
after NFW) density profile with a concentration parame¬ 
ter c = 10. The smaller of the two halos, representative 
of a typical dwarf galaxy (we use the acronym DW), has 
a circular velocity of V 200 = 35km/s and a virial mass of 
M 200 = 1.4 X 10^° M 0 , both measured at the virial radius 
7^200 = 50kpc . The halo is truncated at 112.5 kpc, resulting 
in the total enclosed mass of 2 x 10^° M©. This is essentially 
the same set of parameters we used in the hydrodynamic 
simulations of Dubois & Teyssier (2010) and Teyssier et al. 
(2013). 

The larger halo is chosen to be a typical Milky Way 
galaxy (we use the acronym MW), where we increased 
the circular velocity to V 200 = 160km/s, corresponding 
to a virial radius of i^ 2 oo = 230 kpc and a virial mass of 
M 200 = 1.3 X 10^^ M©. It is truncated at 514 kpc, resulting 
in the total enclosed mass of 2 x 10 ^^ Mq. 

In all other aspects, both initial configurations fol¬ 
low the same prescription as in Teyssier et al. (2013). We 
consider a gas fraction equal to the universal mean value 
/gas = 15%, and the gas density is also following a NFW 
profile. The gas temperature is initialised by solving the hy¬ 
drostatic equilibrium equation. The gaseous halo is set in 
slow rotation around the z-axis, using the angular momen¬ 
tum profile from cosmological simulations and a spin pa¬ 
rameter A = 0.04. The dark matter halo is sampled by 10® 
dark matter particles, whose initial positions and velocities 
were computed with the density-potential pair approach of 
Kazantzidis et al. (2004) and Read et al. (2006). The stabil¬ 
ity of the resulting gas-dark matter equilibrium was shown in 
Teyssier et al. (2013) to be sufficiently good for our present 
purpose. 


3.2 Initial magnetic field 

A fundamental ingredient in any MHD simulation is the 
adopted initial magnetic field configuration. The simplest 
possible choice would be a constant field parallel to one di¬ 
rection, e.g. a vertical uniform field 



However, we argue here that this choice is not appropriate 
for initially concentrated mass distributions in general and 
for cosmological halos in particular. This simple choice re¬ 
sults indeed in a completely uniform magnetic energy distri¬ 
bution. As collapse proceeds, because of the initially peaked 
density distribution, less and less mass is added to the cen¬ 
tral galaxy, especially at late time. Magnetic energy, how¬ 
ever, is still being accreted efficiently, even at late time, and 
artificially added to the central object. We believe it is more 
appropriate, in the context of ideal MHD (frozen-in mag¬ 
netic flux), to consider that the initial magnetic energy fol¬ 
lows closely the initial density distribution, in which case 


the magnitude of the field would scale roughly as 

(14) 

This requirement translates into a more complex field topol¬ 
ogy, and we need to work harder to initialise the field, com¬ 
pared to the uniform, vertical field case. Another important 
property of galactic dynamos is the field parity with respect 
to the system’s mid-plane. To explore possible effects re¬ 
lated to the direction of the vector field, with odd or even 
parity across the mid plane, we consider initially two typi¬ 
cal topologies: dipole-like or quadrupole-like. The dipole-like 
field (we use the acronym D) is defined using the following 
vector potential 


Ad = Bq 


p{r,i 


Po 


2/3 


ref) 


(15) 


where p(r, z) is the initial gas density given by the NFW 
profile and Cf) is the unit vector along the toroidal direction. 
Note that here, coordinates are given in a cylindrical coor¬ 
dinate system centered around the halo center and aligned 
with the rotation axis. The corresponding magnetic field has 
a vertical component which is symmetric with respect to the 
mid plane, while its radial component is antisymmetric. The 
quadrupole field (we use the acronym Q) is defined using the 
vector potential 


An = Bo 


p{r, z) 


1 2/3 


Po 


ZCf) 


(16) 


and has the reversed symmetries. Its vertical component is 
antisymmetric and its radial component is symmetric with 
respect to the mid plane. Figure 1 illustrates the shape of 
the two vector fields. Note that in both cases, the initial 
field strength follows a profile peaked around the center with 
magnitude roughly proportional to as expected. The 

magnetic field is then initialised on each cell face of the AMR 
grid, using a finite-difference approximation of the curl: 


-^initial — V X A 


(17) 


This ensures that the divergence of the magnetic field is 
initially exactly zero (to machine precision), and, thanks to 
the Constrained Transport method, remains zero during the 
entire simulation. 

In the present paper, we would like to explore the purely 
kinematic regime of the magnetic field evolution. This corre¬ 
sponds to very small values of the magnetic field, for which 
there is no back reaction on the flow (the Lorentz force can 
be ignored). We therefore only solve the induction equation, 
which is linear with respect to the magnetic field. The ex¬ 
act value of the parameter Bo is therefore irrelevant, and we 
will always quote magnetic field intensity as a function of 
the initial intensity, or as a function of the average inten¬ 
sity. We will study the saturation regime, and how the field 
reaches equipartition with the thermal and kinetic energies 
of the gas in a companion paper (hereafter Paper H). Close 
to saturation and equipartition, the exact value of the field 
matters a lot, and in this case, the initial intensity plays a 
very important role. Here, however, only the initial spatial 
distribution and the initial topology of the field are impor¬ 
tant, but not its overall initial normalisation. 
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Figure 1. Streamlines of the magnetic field for the two topologies used in the initial conditions. The left panel shows the dipole with 
mid plane symmetry of the vertical component and mid plane antisymmetry of the radial component, while the right panel shows the 
quadrupole field topology with opposite symmetries. 


3.3 Summary of additional physics parameters 

The feedback mechanism, whose details were explained in 
the previous section, can be switched on and off and in 
case of radiative feedback its effective strength can be con¬ 
trolled by changing the surrounding dust opacity parameter 
The dwarf halo simulations were run without any radia¬ 
tive feedback. The only option was to have supernova feed¬ 
back (simulations dubbed ’SN’) or no feedback at all. In the 
Milky Way case, however, adding to those two options we 
also tested two more set-ups with radiation feedback. The 
values tested were one medium scale dust opacity value of 
Av = 5cm^/g (’K-5’) for 70 K dust and a rather opaque gas 
with = 20cm^/g (’K-20’) corresponding to a dust tem¬ 
perature of 140 K (cite Semenov et al. 2003). The feedback 
parameters for the dwarf as well as the Milky Way simula¬ 
tions are given in Table 1 and Table 2. Both configurations 
have the same star formation e* = 1% and supernova feed¬ 
back rjsN = 10% efficiencies. The temperature floor used to 
prevent the gas from fragmenting below our resolution limit 
is given by 

Tmin^n f (18) 

where the critical temperature T* is a cool 100 K for the 
dwarf halo and warm 2000 K in the Milky Way case. 


4 MAGNETIC FIELDS AMPLIFICATION 
THROUGH FEEDBACK PROCESSES 

We will now present the results of our halo simulations, 
where we studied the influence of various stellar feedback 
parameters. This section is organised as follows: First, we 
present our dwarf galaxy simulations, without feedback, 
then using supernovae feedback. Second, we present the 
Milky-Way-sized galaxy simulations. For the latter case, 
supernovae feedback does not differ strongly from the no¬ 
feedback case, although it introduces slightly more turbu¬ 
lence in the gas. Radiation feedback makes however a big 
difference, and we explore two different dust opacities, re¬ 
sulting into two different scenarios for the galactic outhows. 


Table 1. Halo initial conditions parameters (see text for details). 


parameter 

Dwarf 

Milky-Way 

units 

R200 

50 

230 

kpc 

V200 

35 

160 

km/s 

M200 

1.4 X lO^o 

1.3 X 10^2 

M© 

Ax 

18 

84 

pc 

TTlres 

1.5 X 10^ 

1.5 X 10^ 

M© 

m* 

2.0 X 10^ 

5.9 X 10^ 

M© 

T . 

100 

2000 

K 

n* 

14 

4 

H/cc 

e* 

1 

1 

% 

^SN 

10 

10 

% 

Zini 

0.05 

0.05 

Z0 

met. yield 

10 

10 

% 


Table 2. Initial magnetic field topology and feedback parameters 
(see text for details). 


Name 

Topology 

SN Feedback 

Opacity k [cm^/g] 

DW-D 

Dipole 

No 

0 

DW-D SN 

Dipole 

Yes 

0 

DW-Q 

Quadrupole 

No 

0 

DW-Q SN 

Quadrupole 

Yes 

0 

MW 

Dipole 

No 

0 

MW SN 

Dipole 

Yes 

0 

MW K-5 

Dipole 

Yes 

5 

MW K-20 

Dipole 

Yes 

20 


4.1 Dwarf galaxy 

Ah our simulations begin in a similar way, which is the classi¬ 
cal scenario for these cooling halo set-up. The gas, although 
initially in strict hydrostatic equilibrium, loses thermal en¬ 
ergy through radiative cooling. It thus collapses and a cen- 
trifugally supported disk forms from the inside out, thanks 
to the initial angular momentum prohle. In the dwarf galaxy 
case, the disc is relatively thick at hrst: Atomic cooling sets 
a natural temperature hoor around 10^ K. Low tempera¬ 
ture radiative processes (here mostly hne-structure cooling 
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Figure 2. Magnetic energy evolution in the dwarf galaxy simulations without feedback (left) and with supernova feedback (right), with 
dipole (solid lines) and quadrupole (dashed lines) initial conditions. Values are normalised to the initial total magnetic energy in the 
box. The black curves (left) illustrate the shearing amplification Es with compression values of Ec = IOEq (dipole) and Ec = 10^ Eq 
(quadrupole) and a shearing rate of S = (100Myr)“^ , while the black straight line (right) marks an exponential growth exp(t/T) at a 
rate of r = ITOMyr for comparison. Without feedback, magnetic energy rises a few orders of magnitude during the initial collapse but 
is from then only slightly amplified by shearing over the entire simulation time. With feedback, we observe a fast exponential growth of 
magnetic energy over many orders of magnitude. 


of metals) cools the gas further, leading to the formation of 
a thin disc which quickly fragments into dense gas clumps. 
The gas density in these clumps reaches the star formation 
density threshold and the first stars form. 

4 . 1.1 No feedback case 

In absence of feedback, the disc remains very thin, and the 
gas clumps are long-lived. Although our star formation effi¬ 
ciency was set very low (one percent), most of the gas inside 
the dense clumps is converted into stars, after a few disc 
orbital time. The resulting galaxy is very efficient at trans¬ 
forming most of its baryons into stars, which is at odd with 
observed dwarf galaxies in the universe. Moreover, the re¬ 
sulting circular velocity profile is strongly peaked towards 
the centre, although dwarf galaxies circular velocity profiles 
are usually declining towards the centre. 

The magnetic energy evolution can be seen in Figure 2. 
Our new simulations confirm the earlier finding of Wang & 
Abel (2009) and Dubois & Teyssier (2010): During the early 
magnetic field amplification due to the collapse of the gas 
in the first few Myr, the magnetic energy increases as 
where A is the ratio of the gas density after and before the 
collapse. Note that the magnetic field topology matters a lot 
in this early evolution. Without feedback, the dipole config¬ 
uration leads to magnetic reconnection in the mid plane, as 
expected from the antisymmetry in the poloidal field. In the 
quadrupole field configuration, because of the symmetry of 
the poloidal field, the magnetic energy is not affected by 
field cancellation effects. 

After the collapse, one can see in Figure 2 that the 
magnetic energy still grows, but much more slowly. This can 
be explained from field lines being twisted by the differential 
rotation. In this almost perfectly axisymmetric geometry, 
one can indeed approximate the induction equation as (see 


e.g. Dubois & Teyssier 2010) 

dtBr ^ 0 and ~ rBr^rf^. (19) 

The toroidal field grows therefore only linearly with time, 
while the poloidal field (mostly radial) remains constant. 
This results in quadratic time relation of the magnetic en¬ 
ergy 

i;s = • (1 + (S • tf) ( 20 ) 

with the magnetic energy after collapse Ec, which depends 
on the field topology, and the shearing rate S = rdr^. We 
illustrate the contribution of this model for shearing ampli¬ 
fication to Emag for the no feedback simulations in Figure 2 
(left). 

On the other hand, one can see directly from Equa¬ 
tion 19 that a fast, exponential amplification of the field can 
be obtained only if the radial component grows as fast as 
the tangential component. Eor two dimensional, axisymmet¬ 
ric flows like our smooth rotating disk, this cannot be the 
case, according to the famous Zel’dovich and Cowling anti¬ 
dynamo theorems (Charbonneau & Steiner 2012). A further 
inspection of Eigure 2 reveals several spikes in the magnetic 
energy evolution. These are due to collapsing, rotating gas 
clumps, which trigger short episodes of field amplification. 
As anticipated by Wang & Abel (2009), these vortex modes 
do indeed amplify the field locally, but as soon as the clumps 
dissolve in the large-scale rotating flow, so does their mag¬ 
netic energy. 

These clumps also trigger three-dimensional turbulence 
in the gas, thanks to clump-clump interactions (Agertz et al. 
2009). This could in principle increase the magnitude of the 
radial component of the magnetic field, but the induced ef¬ 
fects remain too weak to affect the large-scale dynamo. Eig¬ 
ure 7 shows the velocity dispersion of the dwarf galaxy in 
the no-feedback case: It barely reaches 10% of the tangential 
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Figure 3. Magnetic field components after 1.9 Gyr, as a function of the vertical height relative to the galactic mid plane, normalised 
to BrmSj fhe root mean square (rms) value of the field amplitude. Each component has been computed for each z bin as the volume- 
weighted average value inside each slice (the bin size is 10 pc). Without feedback, the initial mid plane symmetry of the radial component 
is imprinted on the toroidal field, whose peak value is rather high (in units of the rms field). With feedback, all 3 field components are 
equally strong and their typical values comparable in strength to the rms. 


velocity. As a consequence, the corresponding magnetic field 
remains mostly toroidal, as shown in Figure 3. One can also 
clearly see in this Figure that the initial field parity (odd for 
the dipole and even for the quadrupole) has been conserved 
during the collapse and the subsequent shear amplification, 
providing a direct dependence of the final field parity on the 
initial halo field topology. 

4-1.2 Supernovae feedback case 

We now describe our results for the dwarf galaxy with super¬ 
novae feedback enabled. The evolution is drastically differ¬ 
ent, with violent outflows terminating quickly the life of the 
dense star-forming clouds. The resulting star formation rate 
is reduced by one order of magnitude, compared to the no¬ 
feedback case. As shown in Teyssier et al. (2013), the galactic 
circular velocity is now in much better agreement with ob¬ 
served dwarf galaxies, exhibiting a kpc-sized core in the dark 
matter distribution. Star formation also proceeds in succes¬ 
sive starbursts, leading to the ejection of massive quantities 
of gas into a galactic fountain. The gas falls back after a 
dynamical time, triggering a new star formation event. 

The corresponding magnetic energy evolution can be 
seen in in Figure 2. We observe, for both dipole and 
quadrupole initial conditions, a very fast, exponential 
growth with e-folding time of around 200 Myr. The mea¬ 
sured growth rate is therefore quite fast, comparable to the 
rotation rate F ^ Q. Note that the magnetic energy has 
been amplified by almost 18 orders of magnitude, which 
correspond to 9 orders of magnitude in the magnetic field 
itself. Figure 7 compares the rotational velocity V to the 
vertical velocity dispersion a. We have in the feedback case 
l//cr ~ 1, a clear sign of strong turbulence. As a result, field 
lines are violently twisted in random directions, allowing the 
radial and vertical components of the field to reach a similar 
strength than the toroidal component. Figure 3 illustrates 
this isotropy in the magnetic field. The field is now highly 


turbulent, with the largest fluctuation seen on the smallest 
scales (a few cell in size). 

To study further the interplay between the supernovae- 
driven turbulence and the growth of the field, we have com¬ 
puted the power spectra of both the gas kinetic energy and 
the gas magnetic energy in Figure 8. For the former, we 
see the almost immediate onset of a power law power spec¬ 
trum like the one predicted by Kolmogorov (P ~ 
quite typical of supersonic turbulence. The kinetic energy 
power spectrum appears very stable throughout the evolu¬ 
tion, maintained at this high level by supernovae explosions 
and rotational energy. Note that the injection scale of the 
kinetic energy is very large here: it is the size of the en¬ 
tire galaxy. To get an idea of the nature of the forcing of 
the turbulence in this feedback-dominated galaxy, we have 
plotted in Figure 5 a rendering of the gas velocity field. It 
shows large-scale upward and downward motions, together 
with a clear overall rotation pattern. The largest scale at 
which kinetic energy is injected turns out the be roughly 
the halo scale radius r^, for which marks the transition be¬ 
tween the two power law regimes in the dark matter dis¬ 
tribution (from r~^ deep inside to r~^ in the outskirts). In 
what follows, this radius will be identified to the kinetic en¬ 
ergy injection scale, also noted L. Note that in our spectral 
analysis, the global rotation was not removed from the ve¬ 
locity field before computing its Fourier transform, because 
turbulence is clearly dominating. As an additional caveat, 
we also note that both the gas density and the magnetic 
field are far from being homogeneously distributed in the 
box where the spectrum is computed, so that the signals are 
only approximately isotropic and far from periodic. 

The magnetic energy power spectrum, on the other 
hand, is plotted in Figure 8 (bottom). We see here again 
that its amplitude is exponentially growing, while its shape 
remains roughly the same, with P ~ k^^‘^ on the large-scale 
end. The power spectrum reaches a maximum at a scale 
corresponding to 5 cells, then slowly declines as its gets 
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Figure 4. Face-on (top row) and side-on (bottom row) of mass-weighted density projections in the dwarf galaxy simulation at 1.9 Gyr. 
When there is no feedback (left column), the gas builds up a thin and rotationally suported disk. With feedback (right column), we form 
a thick, turbulence-supported disk with strong outflows. 


no feedback 



X[kpc] 



Y[kpc] 


30 

24 

18 

12 

6 

0 

-6 


u 

o 

-12 > 
-18g 
-24^ 
-30 


30 

24 ^ 
18 ^ 
12 ^ 
6 > 
0 G 

05 

-12 > 
-18g 
-24^ 
-30 



Y[kpc] 


Figure 5. Mass-weighted average velocity held maps corresponding to the density projections of Figure 4. Arrows show direction (but 
not strength) of the mean velocity held perpendicular to the line of sight. Colours indicate the strength and direction of the line-of-sight 
velocity component, where blue means approaching and red means receding from the observer. Without feedback (left), gas motion is 
dominated by global rotation (bottom) and the vertical velocity component is weak (up). Feedback (right) drives winds which can be 
seen as a strong vertical component (upper panel) and signiflcant deviations from pure rotation (lower panel). 
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Figure 6. Mass-weighted magnetic pressure maps of the dipole simulation at 1.9 Gyr, normalized to the initial average magnetic pressure. 
Overlaid in black are streamlines of the mean field perpendicular to the line of sight. Without feedback (left), the field is ordered and 
dominated by a large-scale dipole structure (bottom) and a dominant toroidal field in the disk (top). With feedback (right), the magnetic 
field is characterised by random turbulent motions, without preferred direction or large-scale pattern (both right panels). 



0 1 2 3 4 5 6 

radius [kpc] 

Figure 7. Velocity profiles of circular velocity and vertical ve¬ 
locity dispersion in our dwarf galaxy simulations without (top) 
and with feedback (bottom) at 1.9 Gyr. The solid lines show the 
average rotational velocity Vq and the dashed lines shows the ver¬ 
tical velocity dispersion cr^, both as a function of the cylindrical 
radius. 

to the Nyquist frequency of the grid. This is exactly what 
is predicted from Kazantsev’s theory of turbulent dynamos 
(Kazantsev 1968), and confirmed in the forced-turbulence, 
periodic box MHD simulations of Haugen et al. (2004). In 
the present dwarf galaxy simulations, we also obtain a small- 
scale magnetic dynamo, for which the forcing scale would be 


the size of the entire galaxy L ~ 10 kpc, and for which the 
magnetic dissipation scale would be set by the adopted nu¬ 
merical resolution. 

In the small-scale dynamo theory, a critical ingredient 
is the magnetic Reynolds number Rm — VLlr\ which encode 
the magnitude of the small-scale magnetic dissipation. In our 
notation, the velocity dispersion at the forcing scale L is V 
and rj is the magnetic dissipation coefficient. As discussed 
in Brandenburg & Subramanian (2005), exponential growth 
of the field is obtained if Rm > Rcrit, where the critical 
magnetic Reynolds number was observed to be around BO¬ 
SS. In our case, where no explicit magnetic dissipation has 
been included in the induction equation, this translates into 
a critical spatial resolution, beyond which we expect to see 
exponential amplification of the field. 

For this reason, we repeat the dwarf simulations at dif¬ 
ferent resolutions to study their effect. We show in Figure 9 
the time evolution of the total magnetic energy in our dwarf 
galaxy with 3 different maximum resolutions: Ax =36, 18 
and 9 pc. We see a tendency for stronger amplification from 
compression with increasing resolution in all cases because 
the gas can build structures with higher density. Yet, with¬ 
out supernova feedback the subsequent shearing amplifica¬ 
tion does not show a significant impact of resolution in the 
resulting magnetic energy at the end. The feedback-driven 
dynamo, on the other hand, shows a strong growth rate de¬ 
pendence on the maximum resolution. In the low-resolution 
run, we obtain a rather slow amplification, with F ~ 0.4Q. 
This means we are close to or slightly better than the criti¬ 
cal resolution for small-scale dynamo. For our fiducial reso¬ 
lution, we see a fast exponential growth with F ~ Q. In the 
high resolution case, we observe an even faster growth with 


MNRAS 000 , 000-000 (0000) 































A small-scale dynamo in feedback-dominated galaxies 11 




Figure 8. Time evolution of the kinetic and magnetic energy 
spectrum in a central 512^ box and normalised to the total initial 
magnetic energy in the box. The simulation quickly develops a 
k^/^, exponentially growing energy spectrum, typical of small- 
scale Kazantsev’s dynamos, bottlenecked on scales of a few grid 
cells. 

r ~ 2Q. In Figure 10, images of gas density projections are 
plotted to show how the gas structure is resolved at different 
resolutions. With increasing resolution, we see more dense 
substructures and importantly finer details in the flow and 
stronger winds. 

This behaviour can be explained nicely within the 
framework of small-scale dynamos, for which the growth rate 
is determine by the inverse of the eddy turn-over time at 
the dissipation scale namely F ~ v{i)li oc for Kol¬ 

mogorov’s turbulence. At higher resolution, the eddy turn¬ 
over time scale become shorter, so that the growth rate be¬ 
comes larger. In the kinematic phase, well before the satu¬ 
ration phase, when the magnetic energy will reach equipar- 
tition with the kinetic energy, the growth rate of the small- 
scale dynamo is therefore determined by the smallest re¬ 
solved scale of the turbulence, i, and propagates through 
an inverse cascade all the way up to the forcing scale L, 
following the Kazantsev’s power law. 

Since actual microscopic dissipation processes are oc¬ 
curring on very tiny scales, completely unresolved even by 


our highest resolution run, one expects the actual growth 
rate to be even higher than the already fast rates we have 
measured in our numerical experiments. It is therefore rea¬ 
sonable to assume that field saturation will be reached very 
quickly. We however defer the study of the saturation phase 
to a follow-up, companion paper (hereafter Paper II). 

4.2 Milky-Way-like galaxy 

We now describe our results for the Milky-Way-sized halo. 
In this case, the gas quickly cools down to 10^ K, which 
leads to the immediate formation of a very thin disk. Our 
spatial resolution in the Milky-Way case is limited to 84 pc, 
so we can’t allow the gas to cool much below 3000 K, as sum¬ 
marised in Table 1. Nevertheless, the disk is so massive that 
it also quickly fragments into massive clumps that actively 
form stars. 

4-2.1 No feedback and supernovae feedback cases 

In the no-feedback case, the galactic disk develops a massive 
central condensation of stars, which results in a central cir¬ 
cular velocity close to 500 km/s. With supernovae feedback, 
however, we manage to reduce significantly the central bulge 
mass, with a maximum rotational velocity approaching only 
250 km/s (see Figure 15). Note that even with supernovae 
feedback included, we do not reduce the overall star forma¬ 
tion efficiency, and after several Gyr, most of the baryons 
have been converted into stars. This can be explained by the 
relatively low specific energy released by supernovae, mak¬ 
ing it very difficult for the gas to reach the escape velocity of 
the halo (see for example Roskar et al. 2014, for a complete 
discussion). A small galactic fountain sets in, so that the gas 
remains mostly bound to a weakly turbulent disc, as can be 
seen on Figure 15 in the vertical velocity dispersion profile. 

Figure 11 shows the magnetic energy evolution for these 
two cases, and they exhibit the same features as the no¬ 
feedback dwarf galaxy case: early magnetic field amplifica¬ 
tion due to gravitational collapse of the cooling halo gas, 
followed by a weak shear amplification, with some fluctu¬ 
ations associated with fragmenting-clumps-induced vortex 
modes. The corresponding magnetic field topology appears 
as very well organised on large scales, with a dominating 
toroidal component. 

4-2.2 Radiative Feedback case 

Following the methodology explored for the first time in 
Roskar et al. (2014), we now consider a feedback model 
based on radiation from young stars efficiently absorbed by 
dust, and converted into kinetic energy through the infrared 
radiation force. Bear in mind that this model has been de¬ 
signed to be very optimistic, in order to maximise the effect 
of the radiation pressure on dust grains. Realistic modeling 
of the physical underlying processes is not our main objec¬ 
tive here and we rather want to obtain an efficient feed¬ 
back mechanism, launching a strong enough galactic wind, 
and analyse the possible effect of the resulting fountain flow 
on a magnetic dynamo. Images of gas density projections 
are shown in Figure 12. We have modelled our Milky-Way¬ 
like galaxy using two different values for the dust opacity: 
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Figure 9. Time evolution of total magnetic energy in the simulation box for the dwarf galaxy simulation with different maximum 
resolutions without feedback (left) and with supernova feedback (right) and normalised to their initial values. The initial magnetic field 
topology used in all cases is the dipole type. We explore the effect of half and double resolution as compared to the original resolution of 
Ax = 18 pc. The black curve in the left plot illustrates the same shearing amplification as in Figure 2. The black straight lines on the 
right-hand side show exponential growth exp(t/T) with values of r = 500, 150, and 70 Myr respectively. 
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Figure 10. Side-on views of the mass-weighted average density in the central region of the dwarf galaxy simulation with the different 
spatial resolutions as in Figure 9, with a panel size of 20 kpc (top row) and 5 kpc (bottom row). The maximum resolution is increased 
from left to right, starting from Ax = 36 pc (left) to Ax = 18 pc (middle) Ax = 9 pc (right). Increasing the resolution allows one to 
reveal more substructures, finer details in the flow and stronger winds. 


= 5 cm^/g and = 20 cm^/g, which span a range of 
realistic dust temperatures. The latter model gives rise to 
the strongest galactic fountain, and resembles in many as¬ 
pects to the dwarf galaxy case with supernovae feedback. 
The former, lower opacity case appears as less energetic, 
with a weaker wind and slightly smaller turbulence. These 
differences can be seen in the gas images in Figure 13, and 


are expressed quantitatively using the vertical velocity dis¬ 
persion and the rotational velocity profiles (see Figure 15), 
the high opacity, more extreme case giving rise to a quasi- 
spheroidal galaxy, with V/a ^ 2. 

When looking at the corresponding magnetic energy 
evolution on Figure 11, it is interesting to notice that the 
higher velocity dispersion corresponds to the faster growth 
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Figure 11. Magnetic energy time evolution in the Milky Way 
simulations, normalised to their initial values. The dashed line 
marks an exponential growth exp (Tt) at a rate of F = 5 Gyr“^ for 
comparison. Without feedback or only supernova feedback (k = 
0), the energy does not grow after the initial collapse. Increasing 
the effective radiative feedback strength k boosts the growth rate 

r. 

rate. Although both galaxies show an exponential growth 
of the magnetic energy, only the high opacity simulation 
reaches a growth rate as high as F ~ Q. The low opacity case 
only reaches F ~ Q/2. This is again in line with the theory 
of small-scale magnetic dynamos, for which the growth rate 
F ^ v(i)Ii, is proportional to the amplitude of the velocity 
fluctuations. Here again, our turbulent forcing scale is the 
size of the entire galaxy, and the dynamo growth rate is set 
by the the dissipation scale, which in our case corresponds 
to the numerical spatial resolution. 


5 DISCUSSION 

We have performed MHD simulations of feedback- 
dominated galaxies, both dwarf galaxies and Milky-Way¬ 
sized galaxies. We have shown that, if feedback processes 
are strong enough, a small-scale dynamo sets in, with a typ¬ 
ical kinetic energy injection scale L corresponding to the size 
of the entire galaxy (which is in our case is close to the halo 
scale radius r^), and a typical magnetic dissipation scale i 
corresponding in our case to the adopted spatial resolution. 
We have observed an exponential increase of the magnetic 
energy, with growth rate F ^ Q, higher than the galaxy 
rotation rate, possibly much higher if one considers realis¬ 
tic microscopic diffusion processes instead of only numerical 
diffusion. 

Three important aspects are missing in order to apply 
our findings to the origin of galactic magnetic fields: 1- We 
have considered rather idealised simulations of galaxies in 
isolation. 2- We have only described the kinematic phase, 
deferring the discussion of the saturation to a companion pa¬ 
per (Paper II). 3- We have considered feedback-dominated 
galaxies, which are relevant for the high-redshift universe. 
What will happen after this feedback-dominated epoch, for 


quiescent, razor-thin galactic discs ? In this section, we spec¬ 
ulate on possible cosmological consequences of our results on 
the nature of the magnetic field in high redshift galaxies, as 
well as the magnetic strength and topology in lower redshift 
galaxies. 

5.1 Cosmological implications for high redshift 
galaxies 

Although our numerical simulations were not performed in 
a realistic cosmological context, we can still draw conclu¬ 
sions for the cosmic evolution of magnetic fields, assuming 
that the universe is made of a collection of halos of vari¬ 
ous masses, and generalised our results using simple analyt¬ 
ical estimates. For this purpose, we will assume that high- 
redshift galaxies are all dominated by efficient feedback pro¬ 
cesses, so that galactic winds can drive a powerful fountain, 
resulting in a gas rich, turbulence-dominated corona, with 
a size equal to r^, the halo scale radius, that sets the tur¬ 
bulence injection scale. As shown in the previous sections, a 
very efficient small-scale dynamo is likely to develop, with a 
growth rate larger (possibly much larger) than the rotation 
rate of the galaxy. Although a detailed study of the satura¬ 
tion phase is required to study how fast the field will reach 
equipartition (and at what scales, see Paper II), we postulate 
here that each halo reaches equipartition between magnetic 
and turbulent kinetic energy almost instantaneously, within 
a volume set by the halo scale radius Vg. This leads to the 
equipartition value for the field: 

^ = /3gasO-?urb (21) 

For the gas density, we assume that its average value can 
be approximate by the baryon fraction ft = Q^h/^m of the 
total mass halo scale density 

Pgas — fbPs (22) 

and for the turbulence velocity dispersion, we adopt the 
value inferred by our simulations, namely the halo maxi¬ 
mum circular velocity (at Vg) 

(rLb ^ Gax ^ 0.193 X 47rGp«r? (23) 

Using standard redshift-dependant functions for these two 
halo parameters Vg and ps, assuming an average halo con¬ 
centration parameter c = 8 and a baryon fraction fb 0.18, 
we obtain a prediction for the saturated field 

+ (24) 

The resulting magnetic pressure scales as a function of mass 
and redshift exactly like a “viral pressure” in the halo. Note 
that this equipartition field increases quite fast with increas¬ 
ing redshift, but quite slowly with increasing halo mass. 
For a feedback-dominated. Milky-Way-sized galaxy, at red¬ 
shift z = 2, one still predicts a rather strong magnetic field, 
slightly above lOOpG. Although very interesting in setting 
the foundations for a theory of small-scale dynamo in the 
cosmological context, our present discussion remains specu¬ 
lative, in the sense that we do not properly model cosmolog¬ 
ical infah and the associated hierarchical merging of smaller 
structures. This could lead to a dilution of the dynamo- 
amplihed held and lower the growth rate. It is therefore of 


MNRAS 000 , 000-000 (0000) 




14 M. Rieder, R. Teyssier 



-20 -10 0 10 20 
X[kpc] 



X[kpc] 



10^ y 

10^ E 

10 ° ^ 
to 

10 o) 

10 ^ ^ 



Y[kpc] 



Y [kpc] 



U 

U 


X 


to 

c 

O) 

■D 

to 

ru 

U) 


Figure 12. Face-on (top row) and side-on (bottom row) views of mass-weighted average density in the Milky Way simulations at 2.3 
Gyr. Radiative feedback efficiencies are low (left), medium (middle) and high (right). With increased feedback strength, the disk becomes 
thicker and winds stronger. 
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Figure 13. Maps of the mass-weighted average velocity field corresponding to the density projections of Figure 12. Arrows show direction 
(but not strength) of the mean velocity field perpendicular to the line of sight. Colours indicate the line-of-sight component strength and 
direction where blue means approaching and red means receding from the observer. While the velocity is dominated by global rotation 
when feedback is weak (left, SN feedback only), the gas motion develops additional small-scale motion and flows vertical to the disk 
(middle) until it becomes almost isotropic turbulence (right). 


MNRAS 000 , 000-000 (0000) 




















































A small-scale dynamo in feedback-dominated galaxies 


15 





loi 

loj 

10 ^ 

loj 

lOj 

10 ^ 

10 ' 

lOs 

10 ® 

10 ® 

10 ® 

103 

10 ® 


loi 

10' 

10' 

10' 

10' 

10' 

10' 

lOs 

10 ? 

lOe 

10 ® 

10 ® 

103 

10 ® 


I 



a 


Figure 14. Maps of the mass-weighted magnetic pressure corresponding to the density projections of Figure 12 of Milky Way simula¬ 
tions with initial dipole field, normalised to the initial average magnetic pressure. Overlaid in black are streamlines of the mean field 
perpendicular to the line of sight. With weak feedback (left), the magnetic field retains the initial large-scale dipole structure outside the 
disk and a toroidal field inside the disk. Increasing feedback (mid and right) makes the field more small-scale. 



radius [kpc] 


Figure 15. Velocity profiles of circular velocity and vertical ve¬ 
locity dispersion analogous to the one in Figure 7, but for the 
Milky Way simulations with increasing values of the dust opac¬ 
ity K from top to bottom. While the rotational velocity Vq (solid 
lines) becomes smaller, the vertical velocity dispersion cr^ (dashed 
lines) increases. 


primary importance to simulate such feedback-dominated 
galaxy formation models with both MHD and a realistic 
cosmological environment. 


5.2 Transition to quiescent, low redshift galaxies 

The cosmological applications we have derived from a simple 
analytical extension of our numerical results was within the 
context of feedback-dominated, high-redshift galaxies. Low 
redshift galaxies are however quite different. We see in our 
nearby universe many grand design, quiescent disk galaxies, 
with very thin, low velocity dispersion disks, and for larger 
halo masses, we even see red and dead elliptical galaxies. The 
present day universe is therefore very quiet, and strong feed¬ 
back effects are absent, except may be in some very intense 
starbursts, usually triggered by (very rare) merger events. 
Our present methodology, based on isolated, gas rich, cool¬ 
ing halos, does not allow to explore the low redshift universe 
self-consistently, unless one artificially switches off feedback 
processes. This is the strategy we adopt in this section, in 
order to explore the consequence of evolving our simulated 
objects from a feedback-dominated state, to a more quies¬ 
cent state, without strong galactic fountains, resulting in 
much thinner, rotationally supported disks. 

It is indeed very important to estimate how the mag¬ 
netic field, amplified first through a turbulent-driven small- 
scale dynamo, could evolve into a large-scale field, mostly 
driven by rotational shear. The obvious question one might 
ask is: Does the magnetic energy disappear, as the small- 
scale turbulent field reconnects on small scales? How intense 
will the surviving large-scale and mostly toroidal magnetic 
field be, after the galaxy develops into a thin, centrifugally 
supposed disk? For this reason, we decide to re-start our 
dwarf galaxy simulations after 3 Gyr of small-scale dynamo 
amplification, but without any stellar feedback. In the ab¬ 
sence of galactic winds, the turbulent corona rapidly col- 
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Figure 16. Face-on (top row) and side-on (bottom row) views of mass-weighted projections of density (left), velocity fields (middle) 
and magnetic field (right) of suppressed feedback simulation at 3.8 Gyr. After feedback was switched off, the gas collapses into a thin 
clumpy disk. Turbulence becomes weaker, and the magnetic pressure is compressed in the disk with large-scale field lines. 
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Figure 17. Time evolution of the magnetic energy inside the sim¬ 
ulation box of dwarf galaxies with feedback (black) and re-runs of 
the same simulations from 3 Gyr, but with feedback switched off 
(blue), with initial dipole (solid) and quadrupole (dashed) mag¬ 
netic field configuration. When feedback becomes suppressed, the 
exponential growth is halted and the total magnetic energy re¬ 
mains roughly constant. 


lapses back to a thin disk, and a mostly toroidal field ap¬ 
pears after a few rotations. Images for gas density, velocity 
field and magnetic field of the galaxy after it has collapsed 
are shown in Figure 16. 

We see in Figure 17 the evolution of the magnetic en¬ 


ergy versus time for this “suppressed feedback” simulation. 
We see that overall, after the thin disk appears, the mag¬ 
netic energy is for the most part conserved. Two competing 
effects are indeed at work here: gravitational collapse that 
amplifies the field from a corona-diluted state to a thin disk¬ 
concentrated state, on one hand, and magnetic losses due to 
reconnection of mostly small-scale field lines in the mid plane 
of the disk. What is interesting and highly non-trivial, is that 
these two effects basically cancel each other, and that the 
final magnetic energy in the thin disk is the same than the 
initial magnetic energy in the large, turbulence-supported 
corona. 

Large-scale, collapse-amplified magnetic fields have 
replaced small-scale, reconnection-suppressed fluctuations. 
Large-scale modes arise in the small-scale dynamo picture 
because the field is amplified on all scales at the same rate, 
up to the turbulence forcing scale, a well-known property 
of small-scale dynamos. Enough magnetic energy has been 
stored on large scales, so it can survive and compensate for 
reconnection and field cancellation effects during the col¬ 
lapse. 

We would like to stress here again that this new type of 
simulations with suppressed feedback was performed in the 
weak field, pure kinematic regime as well. Since we believe 
that after the feedback-driven, small-scale dynamo phase, 
the field has probably quickly reached equipartition, we need 
to study this transition from thick-corona to thin-disk using 
a properly saturated field, which is the purpose of paper 11. 

Nevertheless, we can also study the topology of the field 
after the disk has settled in a thin, rotation-supported state. 
Figure 19 shows the toroidal, radial and vertical field in the 
disk 2 Gyr after the feedback has been suppressed. The field 
is mostly tangential, with however stronger radial and ver- 
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Figure 18. Velocity profiles of circular velocity and vertical ve¬ 
locity dispersion in the suppressed feedback run, as in Figure 7, 
1 Gyr after the feedback has been switched off. Gas motions go 
back to being rotation-dominated, as in the no-feedback case. 


tical components, compared to the no-feedback case. This is 
in agreement with observations which show pitch angles as 
high as 30° (Patrikeev et al. 2006). The most interesting re¬ 
sults is the topology of the field, which appears quadrupole- 
like, even if we start with only a dipole in the initial con¬ 
ditions. Dipole-like modes present in the large-scale mag¬ 
netic field of the corona have odd parity. They will cancel 
in the mid plane after the gas has collapsed into a thin disk 
(like in the no-feedback case). Quadrupole-like modes, on 
the other hand, have even parity and they will be combined 
non-destructively in the mid plane after the collapse of the 
turbulent corona. One interesting prediction of the scenario 
which we consider here is therefore a systematic quadrupole 
field topology, independently of the initial topology of the 
primordial field. This result is in very good agreement with 
observational data of the Milky Way (Taylor et al. 2009; 
Oppermann et al. 2012), as well as external galaxies (Braun 
et al. 2010; Mao et al. 2012). 


6 CONCLUSIONS 

We have performed MHD simulations of cooling halos, for 
both dwarf and Milky-Way-sized haloes, in the kinematic 
regime where magnetic field are weak enough so that the 
effect of the Lorentz force on the turbulent flow is insignifi¬ 
cant. Using supernovae feedback for dwarf galaxies, as well 
as radiation feedback for large galaxies, we have shown that 
a small-scale dynamo quickly sets in, with the turbulent en¬ 
ergy injection scale roughly equal to the halo scale radius Vg, 
and dissipation scale roughly equal to 4-5 cell sizes. In agree¬ 
ment with small-scale dynamo theory, we observe an expo¬ 
nential amplification of the magnetic energy on all scales, 
up to the injection scale, with a growth rate at least equal 
to twice the disk rotation rate. The growth rate in the kine¬ 
matic phase is set by the adopted numerical resolution. In 
our simulated halo, we need to have a resolution of 20 pc or 
better (for Vg — 3.5 kpc) to obtain a significant growth rate 
(as high as r ~ Q). 

This corresponds roughly to 100 resolution elements 
per turbulent energy injection scale. We believe this reso¬ 
lution corresponds to the critical magnetic Reynolds num¬ 
ber Rm — 35, beyond which magnetic field amplification 
through small-scale dynamo is possible. Note that for weaker 
feedback scenarios, the kinetic energy injection scale is likely 
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Figure 19. Average magnetic field components after feedback 
has been switched off, versus the vertical height relative to the 
galactic mid plane, computed as in Figure 3 and normalized to 
Brms- After collapse, the magnetic field is dominated by its 
toroidal component, but we also see a non-negligible radial com¬ 
ponent. The toroidal component has also developed an even sym¬ 
metry across the mid plane, as in the simulations with initial 
quadrupole and without feedback, even if the initial condition 
was a dipole (top). 


to be smaller, which translates into more stringent resolution 
requirements. This analysis should be kept in mind, when 
one wants to extend this to the cosmological context. One 
difficulty might arise when one considers infall of pristine, 
low magnetic field gas as a source of field dilution. This could 
require a larger small-scale dynamo growth rate for the field 
to be able to increase, and therefore results in even more 
demanding resolution requirements. 

We have shown that a small-scale dynamo, driven by 
strong feedback processes in high-redshift galaxies, could be 
the origin of galactic magnetic fields. This scenario is com¬ 
pletely different than the more traditional large-scale dy¬ 
namo approach, in the sense that the small-scale dynamo 
acts very quickly at amplifying the field up to the turbu¬ 
lent injection scale. The new ingredient is here the fact that 
feedback processes at high redshift are probably energetic 
enough, so that this injection scale is the size of the entire 
galaxy (more precisely the halo scale radius rg). We there¬ 
fore have a small-scale dynamo, together with a large-scale 
forcing, hence enabling the fast amplification of the field all 
the way up to the scale of the entire galaxy. 

We have also shown, using simple numerical experi¬ 
ments with suppressed feedback, that the field can evolve 
into a large-scale toroidal, quadrupole field in a low-redshift 
quiescent state, although it has been amplified by a small- 
scale dynamo. Compared to the large-scale dynamo picture, 
for which the field is amplified in a thin disk over several 
Gyr, the small-scale dynamo picture completely reverses the 
point of view, with a very strong field already in place at 
high redshift (after a feedback-dominated epoch) that slowly 
evolves until the present epoch. This large-scale field could 
even slowly decay, on resistive time scales, and still be large 
enough to account for the observed field strength in nearby 
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galaxies. In this new picture, large-scale dynamos are not 
required anymore to amplify the field from its primordial 
value around to fiG levels. Instead, they are needed 

to maintain the field at a roughly constant level by compen¬ 
sating dissipative losses. 
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